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While theoretical models and simulations of magnetic reconnection often assume symmetry such that the 
magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry 
typically leads to non-ideal flows across the null point. To understand this behavior, we present exact 
expressions for the motion of three-dimensional linear null points. The most general expression shows that 
linear null points move in the direction along which the vector field and its time derivative are antiparallel. 
Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and 
resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point 
motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local 
conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, 
the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the 
null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal 
to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under 
special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. 
The motion of separators cannot be described using solely local parameters, because the identification of a 
particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line. 


I. INTRODUCTION 

Magnetic reconnection^ frequently occurs at and ar¬ 
ound magnetic null points: locations where the mag¬ 
netic field strength equals zero.^“® Magnetospheric null 
points have been identified using multipoint in situ 
measurements as the nulls pass through the spacecraft 
constellation.^ Null points in the solar atmosphere have 
been identified through extrapolation of the photospheric 
magnetic field and morphology in coronal emission.® Nu¬ 
merical simulations of magnetic reconnection and plasma 
turbulence at low guide fields frequently show the forma¬ 
tion and evolution of null points,® as do numerical exper¬ 
iments of typical solar events such as flux emergence.^® 

Two-dimensional, non-degenerate magnetic null points 
are classified as X-type or 0-type depending on the local 
magnetic field structure. If we define M as the Jacobian 
matrix of the magnetic field at the null point, then a null 
point will be X-type if det M < 0, 0-type if det M > 0, 
and degenerate if det M = 0. Magnetic reconnection in 
two dimensions can only occur at null points. 

In three dimensions, the structure of non-degenerate 
magnetic null points is significantly more complex. 
Null lines and null planes are structurally unstable and 
unlikely to exist in real systems.The magnetic field 
structure around a linear three-dimensional null point in¬ 
cludes separatrix surfaces (or fans) of infinitely many field 
lines that originate (or terminate) at the null, and two 
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spine field lines that end (or begin) at the null. A neg¬ 
ative (or type A) null point has separatrix surface field 
lines heading inward toward the null point with spine 
field lines heading outward from the null point. In con¬ 
trast, a positive (or type B) null point has separatrix 
surface field lines heading outward away from the null 
point and spine field lines heading inward toward the 
null point. 

Separators (also known as X-lines by some in the mag¬ 
netospheric community) are magnetic field lines that con¬ 
nect two nulls. Separators that include a spine field line 
are not structurally stable, so separators in real systems 
will almost always be given by the intersection of two sep¬ 
aratrix surfaces. Null points, separatrix surfaces, spines, 
and separators are the topological boundaries that divide 
the magnetic field into distinct domains and are there¬ 
fore preferred locations for magnetic reconnection. 
Three-dimensional magnetic reconnection can also oc¬ 
cur without nulls,especially in regions such as 
quasi-separatrix layers where the magnetic connectivity 
changes quickly. 

Motion of magnetic null points and reconnection re¬ 
gions occurs during any realistic occurrence of magnetic 
reconnection. In Earth’s magnetosphere, X-line retreat 
has been observed in the magnetotail®® and poleward 
of the cusp.®^ At the dayside magnetopause®® and in 
tokamaks,®® the combination of a plasma-pressure gradi¬ 
ent and a guide field leads to diamagnetic drifting of the 
reconnection site that can suppress reconnection. Labo¬ 
ratory experiments frequently show reconnection site mo¬ 
tion and asymmetry, often due to geometry or the Hall 
effect.®'^“®® During solar flares, the reconnection site of- 
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ten rises with time as the flare loops grow and can also 
show transverse motions. 

Theoretical models of magnetic reconnection often as¬ 
sume symmetry such that each magnetic null coincides 
with a flow stagnation point in the reference frame of the 
system. When asymmetry is introduced, there is in gen¬ 
eral a separation between these two points,and in 
some cases a stagnation point might not even exist near 
a null point.In all of these situations, there will gen¬ 
erally be plasma flow across the magnetic null and the 
null will change position. Interestingly, the velocity of 
a null point will generally not equal the plasma flow ve¬ 
locity at the null point.This effect is similar to the 
flow-through mode of reconnection.^^ During asymmet¬ 
ric magnetic reconnection in partially ionized plasmas, 
there may exist neutral flow through the current sheet 
from the weak magnetic field (high neutral pressure) side 
to the strong magnetic held (low neutral pressure) side 
due to the neutral pressure gradient. 

In previous work,^° we derived an exact expression for 
the motion of an X-line when its location is constrained to 
one dimension by symmetry. In resistive magnetohydro¬ 
dynamics (MHD), X-line motion results from a combina¬ 
tion of advection by the bulk plasma flow and resistive 
diffusion of the normal component of the magnetic field. 
In this work, we present exact expressions for the mo¬ 
tion of linear null points in three dimensions and discuss 
the typical properties of the bifurcations of degenerate 
magnetic null points. Section II contains a derivation of 
the motion of linear null points in a vector field. Sec¬ 
tion III uses the results from Section II to describe the 
motion of magnetic null points. Section IV considers the 
local bifurcation properties of magnetic null points and 
provides three examples. Section V contains a summary 
and discussion of this work. 


II. MOTION OF LINEAR NULL POINTS IN AN 
ARBITRARY VECTOR FIELD 


We define x„(t) as the time-dependent position of an 
isolated null point in a vector held B(x,t). We define 
B„(x„(t),t) as the value of the vector field at the null; 
while B„ = 0 for all time, ^ (t) — 0 when the 

null point is moving. We define U to be the velocity of 
this null, 


U = 


dx„ 

d< 


( 1 ) 


The local structure of a non-degenerate null point can 
be found by taking a Taylor expansion and keeping the 
linear terms.The linear structure is then given by 


B = M • (5x, (2) 

where dx = x — x„. The elements of the Jacobian matrix 
M evaluated at the null are given by 


where i is the row index and j is the column index. The 
trace of M equals zero when V • B = 0, and M = (VB)^. 

Next we take the derivative following the motion of the 
null, 

^ + U-VB|^^=0. (4) 

This expression gives the total derivative of the mag¬ 
netic field at the null point using the null’s velocity in 
an arbitrary reference frame. This derivative equals zero 
because the magnetic field at the null by definition does 
not deviate from zero as we are following it. By solving 
for U in Eq. 4, we arrive at the most general expression 
for the velocity of the null point^^“^^ 

U = (5) 

which is valid for vector fields of arbitrary dimension. 
This derivation provides an exact result as long as M is 
non-singular. 

An alternate derivation for Eq. 5 starts from the first 
order Taylor series expansion of B with respect to time 
and space about a magnetic null point, 

^■p> 

B{6x,St) = M • Jx-h ^(5t-kO(||Jxf (g) 

ot 

This first order expansion is valid in the limit of small 
5t and |(5x|. We define (5x„ as the position of the null 
point at 6t. Setting B(Jx„,(5t) = 0 provides a unique 
solution for U = (5x„/(5t, and we again arrive at Eq. 5. 
Unlike the previous paragraph, this derivation uses the 
linearization approximation. Eq. 5 may also be derived 
from the implicit function theorem. 

Equation 5 shows that a null point will move along the 
path for which B and are oppositely directed. The 
null point will move faster if the vector field is changing 
quickly in time or varying slowly in space along this path. 
This exact result for U can be applied to find the veloc¬ 
ity of linear null points in any time-varying vector field 
with continuous first derivatives in time and space about 
the null point. A unique velocity U exists as long as M 
is non-singular. If M is non-singular, then there exists 
exactly one radial path away from the null for which the 
vector field is pointed in a particular direction. 


III. MOTION OF MAGNETIC NULL POINTS 


We next consider the case where B is a magnetic field 
rather than just any vector field. The derivation of Eq. 
5 does not invoke any of Maxwell’s equations. We now 
introduce Faraday’s law. 


aB 

Ih 


= -V X E, 


( 7 ) 


where E is the electric field. By combining Eqs. 5 and 7, 
we arrive at the relation 


My = djBi, 


( 3 ) 


U = (V X E), 


( 8 ) 




where quantities are again evaluated at the null point. 
The first term on the right hand side corresponds to the 
magnetic field being carried with the electron flow ve¬ 
locity, Ve, rather than the bulk plasma flow, the second 
term corresponds to the resistive diffusion of the mag¬ 
netic field at the null, and the third term corresponds to 
the Biermann battery. 


IV. THE APPEARANCE AND DISAPPEARANCE OF 
MAGNETIC NULL POINTS 


FIG. 1. A two-dimensional example showing the motion of 
an X-type null point dominated by resistive diffusion of Bz 
along the z direction. Above and below the null, Bz < 0. The 
negative Bz diffuses along the a direction into the immediate 
vicinity of the null point. At a slightly later time, the mag¬ 
netic held at the current position of the null point will have 
Bz < 0. The negative Bz diffusion cancels out positive Bz to 
the right of the null point, so the resulting null point motion 
is to the right. Reproduced with permission from Ref. 30. 
Copyright 2010 American Institute of Physics. 


which additionally requires continuous first derivatives 
of the electric field in space about the null point. This 
expression does not depend on any particular Ohm’s law, 
and indeed can be applied in situations where there is no 
Ohm’s law. 

Next we consider the resistive MHD Ohm’s law, 

E-bVxB = 7?J, (9) 

where V is the plasma flow velocity and J is the current 
density. The resistivity rj is assumed to be uniform for 
simplicity. Eq. 8 then becomes 

U = V-(10) 

where all quantities on the right hand side are evaluated 
at the magnetic null. This expression requires that B has 
continuous first derivatives in time and continuous sec¬ 
ond derivatives in space about the null point. Null point 
motion in resistive MHD results from a combination of 
advection by the bulk plasma flow and resistive diffusion 
of the magnetic field. Even in the absence of flow, null 
points may still move in resistive situations. The plasma 
flow velocity at the null point does not equal the velocity 
of the null point itself.^® A schematic showing null point 
motion due to resistive diffusion is presented in Fig. 1. 

Equation 8 can also be evaluated using an Ohm’s law 
containing additional terms. For example, we can choose 
our Ohm’s law to be 

E-bViXB = ryJ-b^’^^-^, (II) 
erie erie 

where is the bulk ion velocity, Ug is the electron den¬ 
sity, e is the elementary charge, and pg is a scalar electron 
pressure. For J = eug (Vi — Vg), Eq. 8 becomes 

VUg X Vpe\ 


We next consider the emergence and disappearance of 
magnetic null points, with an emphasis on the instanta¬ 
neous velocity of separation or convergence of the bifur¬ 
cating null-null pair. The local approach taken here com¬ 
plements global bifurcation studies.^® Thus far we have 
only considered non-degenerate null points for which the 
local magnetic field can be described by Eq. 2 using only 
the linear terms in the Taylor series expansion. As long as 
M is non-singular at the null, then there exists a unique 
velocity corresponding to the motion of that null point. 
Non-degenerate null points are therefore structurally sta¬ 
ble and cannot disappear unless M becomes singular.®® 

In contrast, degenerate null points are structurally un¬ 
stable and generally exist instantaneously as a transition 
between different topological states.'^’®® Null points must 
appear or disappear in oppositely signed pairs during a 
bifurcation because the overall topological degree of the 
region cannot change unless a null point enters or leaves 
the domain across a boundary.®® In most situations of 
physical interest, degenerate three-dimensional magnetic 
null points will have rank M = 2 and nullity M = 1 (e.g.. 
Ref. 4). The null space (or kernel) of M will then be 
one-dimensional and corresponds to the eigenvector of 
M with eigenvalue zero. The three eigenvalues must sum 
to zero because of the divergence constraint,® which im¬ 
plies that the two non-zero eigenvalues must either be 
both real and of opposite sign, or both complex and of 
opposite sign.® 

Although the linear representation in Eq. 2 can de¬ 
scribe the magnetic structure surrounding a degenerate 
null point (e.g.. Ref. 6) the region around a bifurcat¬ 
ing null-null pair requires higher-order terms. Third- 
order terms need to be considered only when the first- 
and second-order derivatives both vanish at the null, so 
usually a second-order expansion will suffice. The Tay¬ 
lor series expansion of the magnetic field about a three- 
dimensional null point to second-order in time and space 
is 


B((5x, 6t) 


M • -I- 


dBr. 

dt 


6t - 


1 

2 


6x^ Ha; (5x 
^x^ Hj, (5x 
(5x^ Hz i5x 


-f 


6t 

Y 


(<5x • 

(Sx ■ V)dtBy 

(Sx ■ V)dtB^ 




a^B^ 


0(||<5x||®,Jt®), 


-1 y72 


-1 


2 dt^ 
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where the Jacobian matrix, Hessian matrices, and deriva¬ 
tives are evaluated at the null point. The elements of 
the Hessian matrices are given by Hk,ij = didjBk for 
i,j, k G {a;, y, z}. If the magnetic field is locally continu¬ 
ous, then the partial derivative operators will be commu¬ 
tative and the Hessian matrices will be symmetric. When 
B is constant in time and space, then the fourth and fifth 
terms on the right hand side of Eq. 13 vanish. The po¬ 
sitions of the null points for a given 6t may be found by 
setting B((5x„, St) = 0 in Eq. 13 (or the full expression of 
the magnetic field) and then solving for (5x„. 

The instantaneous velocity of convergence or separa¬ 
tion of a bifurcating null-null pair may be infinite, finite, 
or zero (see Section IV A for examples). Suppose that 
there exists a degenerate three-dimensional null point 
with rank M = 2 and that M has three unique eigenvec¬ 
tors. The first-order directional derivatives of each com¬ 
ponent of B are zero along the one-dimensional null space 
of M. Under most realistic circumstances, the second- 
order directional derivative of B along the null space of 
M will be nonzero. Except in special circumstances, the 
component of velocity along the null space of M of the 
bifurcating null-null pair will be instantaneously infinite. 
Next, consider the two-dimensional subspace that is or¬ 
thogonal to the null space of M. The Jacobian of M in 
this subspace at the null point will be invertible. Conse¬ 
quently, there exists a unique finite velocity within this 
subspace for the two-dimensional null point. Thus, in 
general, the instantaneous component of velocity along 
the null space of a bifurcating null-null pair will be infi¬ 
nite while the components of velocity orthogonal to the 
null space will be finite or zero. 

Next we consider separators that may exist and con¬ 
nect a bifurcating null-null pair. Because these bifurca¬ 
tions cannot change the topological degree of the system, 
the null-null pair will include one negative null and one 
positive null.^’^® Define and E+ as the separatrix sur¬ 
faces and 7 _ and 7 + as the spine field lines of the negative 
and positive null points, respectively. The field lines in 
S_ and 7 + approach the null, while the field lines in E_|_ 
and 7 _ recede from the null. In the neighborhood of a 
linear null, the separatrix surface is given by the plane 
spanned by the two eigenvectors associated with eigen¬ 
values that have the same sign for their real part, and the 
two spine field lines are along the remaining eigenvector.® 

Separators that exist in real systems will almost always 
be given by the intersection of two separatrix surfaces. 
Spine-spine separators may exist if 7 _ and 7 _|_ include the 
same field line. Though spine-spine separators may oc¬ 
cur in some symmetric systems, they are not structurally 
stable and thus can generally be ignored.^® As explained 
above, during the bifurcation of a degenerate null point 
a positive and negative null are formed; hence, spine- 
fan separators can never connect a bifurcating null-null 
pair because such separators connect either two positive 
or two negative nulls. Additionally, not all bifurcating 
null-null pairs will be connected by a separator. 

In most realistic situations, there will not exist a 


straight line separator connecting a bifurcating null-null 
pair as one might intuitively expect (see also Refs. 4 and 
5). Typically, there will exist some angle between the 
separatrix surfaces of each of the two null points in the 
time surrounding the bifurcation. Equivalently, each pair 
of eigenvectors associated with the separatrix surface of 
each null will usually be changing in time, and, in gen¬ 
eral, this evolution will be different for each separatrix 
surface. 

A straight line separator may only be created under 
special circumstances, such as when certain symmetries 
are present. For example, a straight line separator will 
occur if the two nulls from the bifurcating null-null pair 
are both improper nulls (not spiral) and they both share 
the same fan eigenvector which is parallel to the direction 
of motion of the bifurcating nulls. 


A. Bifurcation Examples 

Let us consider a prototypical null point bifurcation of 
the form 


B (x, t) 


(a — z)x + by 
cx — {a + z)y 


+ (5B(x,t), 


(14) 


where a, b, and c are arbitrary real constants with a? + 
be ^ 0. We assume that i5B(x, 0) = 0 so that there exists 
a degenerate null point at the origin with rank M = 2 at 
t = 0. The null space of M at the degenerate null point is 
given by z, which is the eigenvector of M corresponding 
to eigenvalue zero. The remaining eigenvectors of the 
degenerate null point are in the x-y plane and given by 


-b 


-b 

a + Va^ + be 

and 02 = 

a — + be 

0 


0 


which correspond to eigenvalues — {a? + be) and a? + be, 
respectively. We only consider time and space close to 
the bifurcation such that \5Bz\ < |a^ -|- be\. These exam¬ 
ples will elucidate many of the properties of null point 
bifurcations discussed earlier in this section. 


1. First bifurcation example 

Suppose that (5B(t) = —zsgn(t)|t|“ in Eq. 14, with 
a > 1. For t > 0, the third component of B reduces to 
z^ — and two null points exist, but when t < 0, there 
are no null points. At t = 0, a single second-order null 
point appears at the origin as the system undergoes a 
saddle-node bifurcation. For t > 0, the two null points 
are at x„ = [0,0, . The null point with z„ > 

0 will be a positive null if -I- 6c > 0 and a negative 
null otherwise. The null points have velocities of x„ = 

[ 0 ,When 1 < a < 2, the velocity of 
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separation diverges to infinity at t = 0. For the critical 
case when a = 2, the null-point velocities are constant: 
x„ = [0,0, ±1]^. When a > 2, the null-point velocities 
asymptotically approach zero at t = 0. Since the null 
point velocities are purely in the z direction and this is 
also the direction of an eigenvector of the fan planes of 
each of the nulls, the resulting separator is a straight line 
along a; = ?/ = 0 for which < 0- 

Finely tuned examples such as this one are unlikely to 
occur in nature, but show that instantaneous velocities 
that are infinite, finite, or zero are mathematically allow¬ 
able during the bifurcation of a degenerate null point. 

2. Second bifurcation example 

Suppose that (5B(t) = Bt in Eq. 14, where By, and 
B;, are real constants with B^ < 0- By setting B(x„, t) = 
0 , we arrive at this expression for the null point positions 
for t>0, 



FIG. 2. Two improper null points (red and blue spheres) 
resulting from the bifurcation of a degenerate null point with 
+ be > Zn > 0. The fan surfaces of the two nulls (denoted 
by salmon field lines for the positive null at > 0 and light 
blue field lines for the negative null at z„ < 0) intersect to 
yield a curved separator field line (green). The spine field 
lines are orange/purple for the positive/negative null. In this 
example, (a, b, c) = (2, —1,3), B = [1, —1, —1]^, and t = 0.2. 


x„(t) 


-aBxt - bByt =F 

+ be + Bzt 

—cBxt + aByt =F i3y\/ 

+ be-\- Bzt 
±\/—Bzt 


(16) 


The instantaneous velocities of the bifurcating null points 
at t = 0 are given by 


lim x„(t) 
(->■ 0 + 


—aBx — bB 


V 


+ be 

—cBx + aBy 
a? + be 
±oo 


(17) 


The velocity in the x-y plane will be finite except under 
the special circumstance when Bx = By = 0 in which 
case the velocity in the x-y plane will be zero. The in¬ 
stantaneous component of velocity along the null space 
of M is infinite. 

The eigenvectors, eigenvalues, and direction normal to 
the fan for the null points resulting from the bifurcation 
are shown in Table I. The eigenvectors associated with 
each null are not functions of time for this example, but 
the eigenvalues are. The structure of the resulting null- 
null pair depends on the value of be. 

When a? be > z'^ >0 (Case 1), all eigenvalues are 
real so the bifurcation results in a positive improper null 
point with Zn > 0 and a negative improper null point 
with Zn < 0. The separatrix surfaces cannot be parallel 
in such a case, so a separator exists between the two nulls 
as the curved intersection of the two separatrix surfaces 
(see Figure 2 for a typical example). 

When -b 6c < 0 (Case 2), each null has two complex 
conjugate eigenvalues so the bifurcation results in a pos¬ 
itive spiral null point with < 0 and a negative spiral 



FIG. 3. Two spiral null points (red and blue spheres) resulting 
from the bifurcation of a degenerate null point with be < 
0. The spine field lines (thick red/dark blue lines for the 
positive/negative nulls) wrap around each other instead of 
intersecting. The fan surfaces (salmon/light blue lines for the 
positive/negative nulls) are parallel to each other and the x- 
y plane so do not intersect. Hence, no separator connects 
this bifurcating null-null pair. In this example, (a, b, c) = 
(1,1, -2), B = [1,1, -2]^, and t = 0.1. 


null point with > 0. The field lines in the separatrix 
surfaces of both nulls are parallel to the x-y plane, so 
the separatrix surfaces do not intersect to yield a sep¬ 
arator. A spine-spine separator can exist under special 
conditions (e.g., when Bx = By = 0), but under generic 
conditions no separator will exist to connect these two 
newly formed null points. Figure 3 shows an example 
where the spine field lines of each null twist around each 
other before approaching the fan of the other null and 
spiraling away. 

A separator will exist as a straight line between the 
two bifurcating null points if and only if Bx = By = 0 
for both Case 1 or Case 2. 
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TABLE I. The eigenvectors, eigenvalues, and direction normal to the fan for the null points in the second bifurcation example. 



Case 1: 

“h he > 0 

Case 2: 

-b 6 c < 0 


Pos. null {z„ > 0) 

Neg. null {zn. < 0) 

Pos. null {Zn < 0) 

Neg. null (z„ > 0) 

Fan eigenvector, eigenvalue 

Z, 2lZn\ 

z, 212 : 71,1 

61 , \Zn\ + yjo? -b be 

61 , - 2 „ - y/a^ + be 

Fan eigenvector, eigenvalue 

62 , -\Zn\ -t -b be 

61 , \z„\ — y/a? -b he 

62 , \zn\ — y/o? +be 

62 , —\z„\ + y/a? -b be 

Spine eigenvector, eigenvalue 

61 , —\zn\ — y/a? + bc 

62 , \zn\ -\- y/-\-he 

Z, 21 Zn 1 

z, 2 |. 27 j,| 

Direction normal to fan 

[—a + \/a^ + be, —b, 0]^ 

[—a — y/a^ + be, —b, 0]^ 

Z 

z 


3. Third bifurcation example 


In contrast to the first two examples, we now consider 
a magnetic field perturbation that is a function of both 
time and space. We define 5B(x, t) = B(x)t in Eq. 14, 
where and By are real constants and B^ (y) = 3?/ — 2 
is a linear function of y. For the particular case shown 
in Figure 4, (a, b, c) = (1,1, —2). A separator is created 
connecting the bifurcating null-null pair because both the 
eigenvalues and eigenvectors of the nulls evolve in time. 
The separators formed in the case between two bifurcat¬ 
ing spiral nulls are typically long and highly spiralled. 
Solving B = 0 using the above values in Eq. 14 gives the 
following as the null point locations as a function of t, 


x„(t) 


(l ± y/2t - 3tyn) yn 

Vn 

±i/2t - 3tyn 


(18) 


where = [(1 -I- 2t) + \/2it‘^ -|- 4< -|- l] /6t. In the limit 
as t —>■ 0, we must consider the quadratic equation satis¬ 
fied by ?/„, 


^tyl - (1 + 2t)y„ - 2t = 0, (19) 

which implies yn = 0 (and, hence, Xn = Zn = 0) as 
t —>■ 0. Differentiating this quadratic with respect to 
t and then taking the limit t —>■ 0 reveals that = 

—2. It can then be shown that in the limit t —>■ 0, the 
instantaneous velocities of the bifurcating null points are 
x„ = [l,-2,oo]’^. 

V. DISCUSSION 

In this paper, we derive an exact expression for the mo¬ 
tion of linear null points in a vector field and apply this 
expression to magnetic null points. Resistive diffusion 
and other effects in the generalized Ohm’s law allow for 
non-ideal flows across magnetic null points. In resistive 
MHD, null point motion results from a combination of 
advection by the bulk plasma flow and resistive diffusion 
of the magnetic field. These results are particularly rel¬ 
evant to studies of null point magnetic reconnection, es¬ 
pecially when asymmetries are present. Analytical mod¬ 
els of asymmetric reconnection must necessarily satisfy 



FIG. 4. As Figure 3, but here the time derivative of the 
z-component of the magnetic field, Bz, is linear in y. The 
addition of a second-order term in space and time means 
that the fan surfaces (salmon/light blue lines for the posi¬ 
tive/negative nulls) tilt towards each other the instant after 
bifurcation creating a separator which connects this bifur¬ 
cating null-null pair. In this example, (a, fo, c) = (1,1,—2), 
t = 0.1, and 5B(x, t) = \t, 0, Sj/t — 2t]^. 


these expressions. Non-ideal flows at null points allow 
the transfer of plasma across topological boundaries. 

Just as we must be careful when describing the mo¬ 
tion of magnetic field lines,we must also be careful 
when describing the motion of magnetic null points. Null 
points are not objects. A null point is not permanently 
affixed to a parcel of plasma except in ideal or certain 
perfectly symmetric cases. Null points cannot be pushed 
directly by plasma pressure gradients or other forces on 
the plasma, but there will generally be indirect coupling 
between the momentum equation and Faraday’s law that 
contributes to null-point motion. The motion of a null 
point is determined intrinsically by local quantities eval¬ 
uated at the null point. However, global dynamics help 
set the local conditions that determine null-point motion. 

In addition to providing insight into the physics of 
non-ideal flows at magnetic null points and constrain¬ 
ing models of asymmetric reconnection, the expressions 
for null-point motion have several practical applications. 
Locating nulls of vector fields in three dimensions is non¬ 
trivial,^^ but if the null-point positions are found for one 
time, then these expressions provide a method for esti¬ 
mating the positions of null points at future times. When 
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there exists a cluster of several null points, these expres¬ 
sions provide a method for identifying which null points 
correspond to each other at different times. A practi¬ 
cal limitation is that these expressions will often require 
evaluating derivatives of noisy or numerical data (cf. Ref. 
42). However, these expressions provide a test of numeri¬ 
cal convergence and can be used to estimate the effective 
numerical resistivity in simulations of null-point reconn¬ 
ection (compare to Ref. 43). 

Linear magnetic null points appear and disappear in 
pairs associated with the bifurcation of a degenerate mag¬ 
netic null point. The null space of M in these degenerate 
nulls will typically be one-dimensional. Second or higher 
order terms in the Taylor series expansion are necessary 
to describe the structure of a degenerate null point and 
the region between a bifurcating null-null pair. Except in 
special circumstances, the instantaneous velocity of con¬ 
vergence or separation of a null-null pair will typically be 
infinite along the null space of M but with finite com¬ 
ponents of velocity in the orthogonal directions. This 
means that null-null pairs that have just appeared or are 
just about to disappear will not lie next to one another, 
but will always have a finite separation no matter how 
small a time step between frames is taken and regardless 
of whether the field is known numerically on a grid or 
analytically everywhere within a domain. 

Just before or after a bifurcation, a straight line separa¬ 
tor connecting the null-null pair will generally not exist. 
Furthermore, a separator, curved or straight, generic or 
non-generic, will not necessarily connect a null-null pair 
just before or after bifurcation if the nulls involved are 
of spiral type and their separatrix surfaces are parallel. 
The structures of second-order nulls and separators that 
exist very near a bifurcation remain important problems 
for future work. 

In resistive MHD, null points must resistively diffuse 
in and out of existence. In the reference frame of the 
moving plasma, a necessary condition for a degenerate 
null point to form is that the resistive term in the in¬ 
duction equation, r^V^B, be antiparallel to the magnetic 
field at the location of the impending degenerate null. 
This places physics-based geometric constraints on when 
and where bifurcations are allowed to happen. 

We may consider whether or not a similar local analysis 
can be performed to describe the motion of separators. 
Consider a separator that connects two magnetic null 
points. Suppose that a segment of this separator exhibits 
non-ideal evolution. Along the remainder of its length, 
the magnetic field in the vicinity of the separator evolves 
ideally. At a slightly later time, the field line that was the 
separator will, in general, not continue to be the separa¬ 
tor between these two null points despite the locally ideal 
evolution. The motion of separators therefore cannot be 
described using solely local parameters. However, it may 
be possible to derive an expression for the motion of a 
separator by taking into account plasma flow and connec¬ 
tivity changes along its entire length. Such an approach 
would provide insight into the structural stability of sepa¬ 


rators and separator bifurcations,^"^ as well as the nature 
of plasma flows across topological boundaries. This lat¬ 
ter aspect is fundamental to the basic physics of three- 
dimensional reconnection; indeed, an early definition"^® 
states that reconnection is “the process whereby plasma 
flows across a surface that separates regions containing 
topologically different magnetic field lines” (see also Ref. 
18). We are investigating the problem of separator mo¬ 
tion in two and three dimensions in ongoing work. 

There exist numerous additional opportunities for fu¬ 
ture work. Our results take a local approach; conse¬ 
quently, numerical simulations are needed to investigate 
the interplay between local and global scales during null- 
point motion and bifurcations. Numerical simulations 
can be used to investigate how null points diffuse in and 
out of existence in non-ideal plasmas and how separa¬ 
tors behave during bifurcations. If the flow field and 
magnetic held are well diagnosed in space or laboratory 
plasmas, the expressions for null-point motion may be 
used to provide constraints on magnetic field dissipa¬ 
tion. Equation 10 offers another opportunity to measure 
the plasma resistivity in the collisional limit. Dedicated 
laboratory experiments offer an opportunity to investi¬ 
gate plasma flow across null points (and other topolog¬ 
ical boundaries) as well as null point and separator bi¬ 
furcations. Finally, many results exist in the literature 
outside of plasma physics on bifurcations of vector fields 
and topology-based visualization of vector fields. While 
communication across disciplines is hindered by differ¬ 
ences in terminology,^® the application of this external 
knowledge to plasma physics will likely lead to improved 
physics-based understanding of these processes. 
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